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Abstract 

We calculate the spectrum of excitations in strong-coupling lattice QCD in a background of 
fixed baryon density, at a substantial fraction of the saturation density. We employ a next-nearest- 
neighbor fermion formulation that possesses the SU(Nf) x SU(Nf) chiral symmetry of the con- 
tinuum theory. We find two types of massless excitations: type I Goldstone bosons with linear 
dispersion relations and type II Goldstone bosons with quadratic dispersion relations. Some of the 
type I bosons originate as type II bosons of the nearest-neighbor theory. Bosons of either type can 
develop anisotropic dispersion relations, depending on the value of Nf and the baryon density. 
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I. INTRODUCTION 



In a previous paper we constructed a framework for calculating the effects of a back- 
ground baryon density in Hamiltonian lattice QCD at strong coupling. We used strong 
pling perturbation theory to write an effective Hamiltonian for color singlet objects 
At lowest order we obtained an antiferromagnetic Hamiltonian that describes meson physics 
with a fixed baryon background distribution. (Baryons move only at higher order.) The 
Hamiltonian was then transformed to the path integral of a nonlinear a model. The latter 
is most easily studied at large N c . 

The global symmetry group of the a model depends on the fermion kernel of the lattice 
QCD Hamiltonian. For Nf flavors of naive fermions we get an interaction between nearest- 
neighbor (NN) sites that is invariant under U(N) with 

N = ANf. (1.1) 

This symmetry is too large and is indicative of species doubling. We add next-nearest- 
neighbor (NNN) interactions to the kernel and reduce the symmetry to 

U(N f ) L x U(N f ) R , (1.2) 

which is almost the symmetry of the continuum theory The unwanted U(1)a is inevitable 
if one starts with a local, chirally symmetric theory |^|. It can easily be broken by hand in 
the a model and we ignore it. 

In 0] we studied the NN theory and found that the ground state breaks U(N) sponta- 
neously. The breakdown pattern depends on the baryon density. In |5] it was shown that 
the excitations of the NN theory divide into two types: type I Goldstone bosons with linear 
dispersion relations and type II Goldstone bosons with quadratic dispersion relations. These 
excitations fit the pattern described by Nielsen and Chadha 0| and studied by Leutwyler [?} . 
Type II Goldstone bosons, typical of ferromagnets, are prominent in work on effective field 
theories for dense QCD [1, H| . 

The NNN interactions may be treated as a perturbation that removes some of the global 
degeneracy of the NN vacuum. In Ref.E3 we found that for all baryon densities studied, the 
ground state breaks the NNN theory's axial symmetries. In all cases with nonzero baryon 
density, the discrete rotational symmetry is broken as well. In this paper we investigate how 
the NNN interactions affect the Goldstone boson spectrum, completing our picture of the 
lattice theory that has the symmetry of continuum QCD. 



II. NON-LINEAR a MODEL 

We give here a brief description of the elements comprising the o model. More details 
may be found in []]]. 

The a field at site n is an N x N hermitian, unitary matrix given by a U(N) rotation of 
the reference matrix A, 

a n = U n AUl (2.1) 

with 
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The a field thus represents an element of the coset space U(N)/[U(m) x U(N — m)]. The 
number m can vary from site to site and is determined by the local baryon number B n 
according to 

m = B n + N/2. (2.3) 

The Euclidean action is 



S 



N 



dr 



]T Tr AC/t d T U n + ^ E Tr 



J> 



^Tr (a n aiVn+2iai) 



(2.4) 



Here is the 4x4 Dirac matrix, times the unit matrix in flavor space. The NN term is 
invariant under the global U(N) transformation U n — > VU n (or a n — > Va n V^) while the 
NNN term is only invariant if V^a^V = oti for all i. This restricts V to the form 



V = e W [i(6 a v + l5 9 a A )\ a } 



(2.5) 



where A a are flavor generators. This is a chiral transformation in U(Nf) x U(Nf). [The 
£7(1) corresponding to baryon number is realized trivially on <j n .] 

The NNN term couples (discrete) spatial rotations to the internal symmetry, viz., 



R^a-n'R, n = Tin. 
Here TZ is a 90° lattice rotation and R represents it according to 



R 



cxp 



.It 1 

I— 




°)1 


. 4 ' 




a i J. 



® 1 N 



(2.6) 
(2.7) 

If the NNN fermion kernel is taken to be a truncated SLAC derivative [llj, then both 
couplings Ji and J 2 ar e positive, and J 2 = Ji/8. If we argue, however, that the strong- 
coupling Hamiltonian is derived by block-spin transformations applied to a short- distance 
Hamiltonian, then we cannot say much about the couplings that appear in it. We will 
assume that couplings in the effective Hamiltonian fall off rapidly with distance; indeed we 
will assume that 

< J 2 < Ji/JV c . (2.8) 

This means that we take as our starting point the (globally degenerate) vacuum determined 
in for the NN theory with 0(1/ N c ) corrections. The NNN interaction is a perturbation 
on this vacuum and its excitations. 



III. THE GROUND STATE 

Here we give a short reprise of the results of 0, for the ground states of the NN 
and NNN a models with a uniform baryon density B n = B > (i.e., a uniform value of 
m > N/2). 

The overall factor of N c in Eq. (|2.4jl allows a systematic treatment in orders of 1/N C . 
In leading order, the ground state is found by minimizing the action, which gives field 
configurations that are time-independent and that minimize the interaction. Minimizing 
the NN interactions results in a locally degenerate ground state: We assign o = A on the 
even sites and let the a field on each of the odd sites wander freely in U(m)/[U(2m — N) x 
U(N — m)], a submanifold of U(N)/[U(m) x U(N —m)]. Since the odd sites are independent, 
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the degeneracy is exponential in the volume. 1 In Ref.Qwe showed that 0(Ji/N c ) fluctuations 
generate a ferromagnetic interaction among the odd sites, causing them to align to a common 
value ( "order from disorder" ) . 

The resulting ground state has a Neel structure. The even sites break U(N) 
to U (m) x U (N — m) and then the odd sites break the symmetry further to 
U(2m — N) x U(N — m) x U(N — m). We can write explicitly 



Ceven UA. even U^ U 

a odd = UA odd U* = U 



1 IjTt 









— ljV-m / 


^ l-2m- 


-N 








-1n- 


-m 


V o 





ljV-m 



{7 f . (3.1) 



The matrix {7 G U(N) represents the global degeneracy due to spontaneous symmetry 
breaking. 

We showed in Ref. [H3 that the NNN interactions partially remove this global degeneracy. 
We made the ansatz 

and showed that it yields a ground state. When m > 3N/4, the matrix u is free to take any 
value within U(N/2), but a U(2m - 3N/2) x U(N - m) X U(N - m) subgroup of U(N/2) 
acts trivially in Eq. ()3.1|) (i.e., matrices in this subgroup give the same field configuration 
as choosing u = 1). Vacua that are associated with different nontrivial choices of u are in 
general inequivalent, and give different realizations of the U(Nf) x U(Nf) symmetry of the 
theory. Since these vacua are not related by symmetry transformations, there is nothing to 
prevent lifting of the degeneracy in higher orders in 1/N C . In the sequel, we set u = ljv/2- 
This gives the vacuum with the largest symmetry accessible via the ansatz (|3.2j) . 

For m < 3N/4, u was found numerically by minimizing the NNN energy (|2.4jl . In view 
of what happens for m > 3N/4 this may be only one point in a degenerate manifold of 
ground states. 2 We emphasize that the degeneracy of these vacua is not related to the 
global U(Nf) x U(Nf) chiral symmetry. It is an accidental global degeneracy of the ground 
state. 

The symmetries of these ground states are summarized in Table HI reproduced from (lp| . 
In general, both chiral symmetry and discrete lattice rotations are broken; in some cases a 
symmetry under rotations around the z axis survives. Note that if we remove the unphysical 
axial i7(l) symmetry from the a model, all its realizations will also drop from TableHJ namely, 
there will be no unbroken axial {7(1) symmetries (third column) and no Goldstone bosons 
corresponding to a broken axial {7(1) (fourth column). 



1 This degeneracy is not removed by the NNN interactions. That is why we consider the 0(1/N C ) corrections 
first. 

2 We have found, in fact, one case where a different ansatz gives a more symmetric ground state than 
Eq. This is the case (N = 12, m = 8), i.e., (N f = 3, B = 2). 
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TABLE I: Breaking of SU(Nf)L x SU(Nf)n x U(1)a for all baryon densities (per site) accessible 
for Nf < 3. 



N f 


\B\ 


Unbroken symmetry 


Broken charges 







— 


1 


1 


1 


— 


1 




2 


U(1)a 










SU{2) V 


4 




1 


u(i) h 


6 


2 


2 


SU(2) V 


4 




3 


u(i) h 


6 




4 


SU(2) L x SU(2) R x U(1) A 










SU(3) V 


9 




1 


U(1) Y x SU{2) V 


13 




2 


U(1) Y 


16 


3 


3 


SU(3)v 


9 




4 


U(l) h x U(1) Y 


15 




5 


U(l) h x U(1) Y x 


14 




6 


5C/(3)i x SU(3) R x 






IV. SPECTRUM OF EXCITATIONS 

The Goldstone bosons of the NN theory were discussed in As mentioned, they divide 
into two types. There are 2(N — m) 2 bosons of type I with u ~ J\\k\ at low momenta; 
these are generalized antiferromagnetic spin waves (and are the only excitations at zero 
density). There are also 2 (2m — N)(N — m) bosons of type II, that derive their energy 
from quantum fluctuations in 0(1/N C ). These are generalized ferromagnetic magnons with 

u~(j x /N c )\k\ 2 . 

The two types of Goldstone bosons belong to different representations of the unbroken 
subgroup U(2m — N) x U(N — m) x U(N — m). This means that they cannot mix to any 
order in 1/N C . The type I-type II classification is robust in the NN theory. 

Now we calculate the effects of the NNN interactions on the spectrum. In view of Eq. (|2.8J) , 
the NN contributions to the propagators, found in 0], remain unchanged. In particular we 
can take over the self-consistent determination of the self-energy of the type II bosons. We 
need consider the NNN contributions to the propagators in tree level only. We proceed to 
calculate these for m > 3N/4. In these cases, the calculations simplify (much as in ^|) and 
we perform them analytically. 3 We believe that the spectra of the other cases have similar 
features. 

In the NN theory the a fields represent fluctuations around the vacuum (|3.1|) with U = 1. 



3 An exception is (N — 12, m= 10), where we have no analytic solution. See below. 
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We parametrize them [5|] as 



(l-2 X X ] -2x7r f 

-27TX f 1 - 27T7T t 

V -2S*t ~2S^ 



-2x5 \ 
-2ttS 
•1 + 20t0 / 



(4.1) 



and 



O'odd 



/ 1 " 2xx f 
V 2vr X t 



-2 X S 2 X 7rt \ 
-1 + 2</>V 257r f 



27T.S 



(4.2) 



1 - 27T7T t / 



Here <p is an m x (AT — m) complex matrix field written as 



X 

7T 



(4.3) 



and 5 = y 1 — 0t</>. The field ir is an (AT — m) x (AT — m) complex matrix, representing the 
type I Goldstone bosons, x is a (2?7i — AT) x (AT — m) complex matrix that represents the 
type II bosons. If = 0, we have o" eve n,odd = A everii0 dd, which is the ground state of the NN 
theory. We adapt Eqs. (j4.1j) - ()4.2j) to the NNN theory by rotating them, 

UaU\ 



a 



(4.4) 



with U as given in Eq. (|3.2j) . Now = corresponds to the ground state of the NNN theory. 

We substitute Eqs. (I4.1j) - (|4.2|) into the action (|2.4|) . The rotation U disappears from the 
kinetic term and from the NN interaction — they are both U(N) invariant. This means that 
the bare spectra found in when U was absent, remain intact. 

We write the NNN energy as 



EL 



N C J 2 



(4.5) 



aNi 



where a = (even, odd) and N denotes a site on the corresponding fee sublattice. We rescale 
0^0/ y/Nc and expand Eq. (|4.5jl to second order, 



E nnn = ^o + ^V^^Tr^A^^A^ + Ag 



4 



) 

N+2i 



aNi 



+T E Tr ^ A ^^ A SU 2 i + yE Tr *Aa* (A« + Ai : 

4 aNi 4 aiVi 

/ 1 \ 

+0 



(2) 

jv+2-; 



(4.6) 



We have defined ai = U^aJJ , and A^ 1 ' 2 - 1 correspond to the linear and quadratic deviations 
of the a fields from their ground state values. The latter are given by 

( 

Ai 2 ) 



A (D 



-2 X \ 
-2tt 

\-2 X ] -2vrt j 



I -2 X X ] -2^t \ 

-2liX ] -27T7T t 

0+2 



(4.7) 
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on the even sites, and 

/ -2 X \ 
= V -2tt 
\ -2 X f -27rt y 

on the odd sites. Here 



Al 2 > = V 



(-2 X X ] -2%7rt \ 

-27TX t -27T7T t 

v +20V y 



V f (4.8) 



f l o o\ 

1/= 1 (4.9) 

v° - 1 °; 

is the matrix that rotates A even to A D dd- It is easy to show that the terms linear in A*- 1 ) 
vanish. 

In view of the block structure of A ev en,odd and of U, as given in Eqs. (|H.lj) - (|H.2|) . and of 
ai, it is convenient to decompose x f° r m > 3N/4 as 



X 



X2 



(4.10) 



Here xi has N/2 rows, and X2 has 2m — 3N/2 rows. Both have N — m columns. Substituting 
into Eq. (J4.6)) and omitting the ground state energy we find that the 0(1) contribution of 
the NNN energy depends only on xi- The tc and X2 fields do not enter the NNN energy at 
this order and remain of type I and of type II, respectively. 
We now define a new N x N matrix x that contains only Xi, 



X 



( xi \ 


\xl o o ) 



(4.11) 



and use it to write 



£ |4Tr [adN^x^i + * (Vx° N V^ a t (Vx° N+2 ^ 



jvefcc 



-4Tr 



a 



AeaiAe (xn? + «iA aiA (yx° N V ] 



(4.12) 



Next we expand x = X V T V , where x v are rea I an d ^ v are the hermitian generators of U(N), 
normalized to 



tr 



8™ . (4.13) 

The form ()4.11|) of x implies that x v 7^ f° r those generators whose elements (J? v ) a p are 
nonzero for a G [1, N/2] and (3 G [m + 1, N] or vice versa. Thus 



unn — 



X^(iV a ) W X^'+Ex^(M ( 



N+2i 



(4.14) 



where 

(iV e ) W 
(iV ) W 



- J 2 Ei Tr 
-J 2 Ei Tr 



VTT^yt^AaOiAa 



(M e 



J 2 Tr 
J 2 Tr 



VTW^aiVT^'V^ai 

(4.15) 
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With a Fourier transform, 



AN 



Ea 



J N 8 ' 

we write the energy in momentum space as 

£nnn = £ X? [ N e + Nj + M e (k) + M\(k 

+X°\ \N + Nj + M {-k) + Ml(-k) 



(4.16) 



feesz 

w>0 



xl 



(4.17) 



HereM a (fc) = Z l M ai e ik *. 

The NN action, including the time derivative and 0(1/N C ) self-energy, was written down 
in |j| in terms of the Fourier transform \ of the (2m— N) x (N—m) matrix field x [Eq. (|4.10|) ]: 



S nn = £ Tr (iw - Si >fe ) XkXl + (-^ - s i,fc) X%X C - 



(4.18) 



The self-energies E a are of order Ji/N c and depend on N and m. We set X2 = and repeat 
the steps leading to Eq. ()4.17|) to write S^n in terms of xl, 



Snn = £ X^KeeXt + X°l^ooX\ + xf KoX"- 



•fc> 



(4.19) 



with the matrices -K" ee , i^ 00 , and K eo given by 



(^e 



'/'/ 



— Tr 

2 

iui . 



-Tr 

2 

S 2 ,feTr 



A^r" 



(4.20) 

(4.21) 
(4.22) 



Equation ()4.19|) is to be added to Eq. (|4.17jl to give the quadratic action of the type II 
Goldstone bosons xi- 

Diagonalizing the quadratic form is straightforward but tedious. As noted above, only 
a subset of the generators T v of U(N) appear in the expansion of Eq. (|4.1ip . Since Xi has 
dimensions N/2 x (N — m), there are N(N — m)/2 pairs of generators in the sum, which we 
write (similar to the Pauli matrices a x , a y ) as T^, T^, with rj = 1, . . . , N(N — m)/2. Their 
coefficients are similarly written as x& Xy- Thus for each k, r\ we have 



xl 



( x e 2 \ 

x e l 

A.yk 
or/* 
Xx,—k 
i on* 
V Xy '-k ) 



and the action is 



s = J2J2xlG- k 1 Xk , 

w>0 k 



(4.23) 



(4.24) 
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The inverse propagator is 



Gi. 1 



/ -Si jfc + n e + m e 

UJ 

V o 



where 



-u -S 2 ,fc \ 

-Ei jfe + n e - m e S 2 ,fc 

-Si >fc + n + m —uj 

S 2ifc u + n -m J 



5 W 



2J 2 $> [r^H' 



cos h. 



(4.25) 

(4.26) 
(4.27) 



Si )2 and uj contain a factor of S^. If we write the 4x4 matrix G k 1 in terms of 2 x 2 blocks, 

/ A B\ 



Gk 1 



C D 



(4.28) 



then its determinant is easily calculated via 

\G^ l \ = \C\ \B-AC~ 1 D\ 

uj 2 + (n e + m e - Si)(n + m - Si) - Ei 
oj{n + m — n e + m e ) 



u{n e + m e — n Q + m ) 
uj 2 + (n e - m e - Si)(n c - m - E x ) - 



(4.29) 



For N = ANf < 12 (and m > 3N/4), the matrices n e ,n ,m e ,m all commute, except 
for the case (N = 12, m = 10). Dropping this last from consideration, we are left with the 
values of (Nf,B) listed in Table Hi) For each case, the simultaneous diagonalization of n 6i0 
and m efi gives the eigenvalues shown. The zeros of the determinant ()4.29j) determine the 
spectrum u>(k), giving (after uj — »■ iuj) 



to 



Si 2 - E 2 2 + - (n 2 - m 2 + n 2 - m 2 ) - Ei(n e + n, 



±< 



Si(n - n e ) + -(n 2 -m 2 e -n 2 + m 2 ) 



i 2 



1/2 



+E 2 2 



[m e + rn ) 2 - (n e - n ) 2 



(4.30) 



Because of the symmetries of Eq. fl4.30j) , the spectra of the various cases shown in Table |H] 
fall into four classes. We examine each class in turn. 

Class 1: Here there is no contribution at all from the NNN interaction. As shown in 
0, Si j2 are proportional to J\/N c and for small momenta they are quadratic in |fc|; the 



same holds for the NN energy y Ex 2 — E 2 2 . These fields remain isotropic type II Goldstone 
bosons as in the NN theory. 
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TABLE II: 


Simultaneous 


eigenvalues of n t 


;0 and m e>0 


(in units 


of 2J 2 ) for all 


values 


of N f 


and B considered. 


The 


resulting spectra fall into four classes. 


Here v z = cos 


k z and 


v± = 


2 sin 


)sin( fa 
















N f 


B 




n 


m e 


m G 


multiplicity 




class 


1 


1 


1 


2 




v±_ 






4 






2 


1 




Vz 






4 


2 


2 


1 


2 


v z 


v± 


x3 




4 






2 


1 


v± 


Vz 


x3 




4 






1 


2 


-v z 


-v± 






4 






2 


1 


-v± 


-V z 






4 


2 


3 





1 





V z 






2 






1 





v z 









2 









2 





v± 






3 






2 





U_L 









3 


3 


3 


1 


2 


^2 


v± 


x6 




4 






2 


1 


U ± 


Vz 


x6 




4 






1 


2 


-V z 


-v± 


x3 




4 






2 


1 


-v± 


-V z 


x3 




4 


3 


5 














x2 




1 









1 











2 






i 
l 














2 









2 





U_L 






3 






2 





v± 









3 



C/ass Fields that correspond to the minus sign in Eq. (|4.H()jl remain type II, but with 
anisotropic dispersion laws of the form 



The plus sign in Eq. ()4.30|) gives a linear dispersion law, again anisotropic, 

lu 2 = 4ci J 2 \k 2 x + k 2 y + (1 + b5)kf 



(4.31) 



(4.32) 



The anisotropy in both cases is proportional to the ratio 5 = J2/(12 Ji/N c ) of NNN to NN 
couplings. The coefficients c and c\ are defined as 



d 



Cl 



Si 



2 -s 2 2 



fc=0 



> 0. 



They are of order Ji/N c . The coefficients 

_ N c c\ - c 2 
G 6Ji 2 Cl c 2 



fc=0 



and 6 



6Ji Ci 



(4.33) 
(4.34) 

(4.35) 
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are of order 10 2 for the cases at hand. 

Classes 3 and 4 : Taking k = in Eq. ()4.30j) we find that the fields that correspond 
to the plus sign get a mass equal to 2J 2 . This is a result of the explicit breaking of the 
U(ANf) symmetry by the NNN interaction terms; these particles are no longer Goldstone 
bosons. The massless bosons in Class 3, corresponding to the minus sign, are type II bosons 
described by 

uj 2 = cjk 4 . (4.36) 

(This is different from Class 1 where uj 2 = c 2 k 4 .) The massless bosons in Class 4 are again 
anisotropic, obeying Eq. (|4.32p . and are of type I. 

These dispersion relations are correct to 0(5) for momenta of 0(5) or smaller. In all cases 
the dispersion relation to 0(1) for k 2 ^> 5 is quadratic and isotropic, unchanged from the NN 
result presented in 0. Since 5 is a small parameter, this means that in most of the Brillouin 
zone the propagator maintains its NN form. This is the reason why the self consistent 
calculations in p| that yield £1.2 do not change when we add the NNN interactions. 



V. SUMMARY AND DISCUSSION 



In this work we have studied the non-linear sigma model derived in 0] for the description 
of lattice QCD with a large density of baryons. The model has NN and NNN interactions. 
Building on the results given in [1 
NNN ground state presented in |l 



and p for the NN theory, and on the study of the 
, we have determined the dispersions relations for the 
Goldstone bosons in the NNN theory. 

We find that the physics of the NN theory is mostly undisturbed by the NNN interaction. 
At leading order, the properties of the type I bosons (tt) and of some of the type II bosons 
(% 2 ) do not change. The type II bosons grouped in the Xi field suffer a variety of fates, 
falling into four classes that appear for different values of Nf and m. Class 1 bosons are 
unaffected by the NNN perturbation. Class 2 bosons split into type I and type II, and all 
gain anisotropic contributions of 0(5) to their energies. Some of the Class 3 bosons become 
massive while other remain unaffected. In Class 4, some become massive while the others 
become anisotropic type I bosons. 

The symmetry of the theory, in all cases, is severely broken by the NNN terms — from 
SU(4Nf) to SU(Nf) x SU(Nf) x U(1)a- Not surprisingly, a simple count shows that the 
total number of massless real fields is far greater than the number of spontaneously broken 
generators of SU(Nf) x SU(Nf) x U(1)a, as shown in Table |1J The particular NNN inter- 
action we use is simply unable to generate masses in lowest order for many of the particles 
unprotected by Goldstone's Theorem. This is partly reflected in the accidental degeneracy 
of the ground state, which we mentioned below Eq. ()3.2j) . Just as this degeneracy should 
be lifted in higher orders in 1/N C [beginning with 0(J2/N C )], the corresponding massless 
excitations should develop masses. The only particles protected from mass generation are 
the minimal number needed to satisfy Goldstone's theorem (or the Nielsen- Chadha variant). 

Another effect that is missing is the mixing of type I and type II Goldstone bosons, which 
is certainly permitted when the NNN interaction is turned on. In (if we proved that such 
a mixing is forbidden in the NN theory, since the two types of boson belong to different 
representations of the unbroken subgroup. The classification in the NNN theory is less 
restrictive, and permits mixing of the bosons. Whether mixing occurs is a dynamical issue 
that can be settled only by calculating to higher order in 1/N C . 
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To conclude we note that other recent work |8|, |9j on the high density regime of QCD — in 
the continuum — also predicts type II Goldstone bosons and anisotropic dispersion. There, 
the starting point is an effective field theory that describes the low energy dynamics of QCD 
with nonzero chemical potential \x. For /i ^ 0, Lorentz invariance is broken, and the field 
equations become nonrelativistic. This leads to the emergence of type II Goldstone bosons. 
In addition, the ground state in |2| can support a nonzero expectation value of vector fields. 
This breaks rotational symmetry and makes some of the dispersion relations anisotropic. 
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